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Abstract 



In a reasonably self-contained and explicit presentation we illustrate 
the efficiency of the Feynman-Kac formula for the rigorous derivation 
1^ ■ of three inequalities of interest in non-relativistic quantum mechanics. 
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"As it was, as soon as I heard Feynman describe his path integral approach 
| to quantum mechanics during a lecture at Cornell everything became clear 

at once; and there remained only matters of rigor. " [1] 
• Jh . 

"Hardly a month passes without someone discovering yet another 
application [of the Feynman-Kac formula] . Remarkable and, of course, 

highly gratifying. " [2] 

Mark Kac (1914-1984) 

1 Feynman-Kac formula 

In order to model the phenomenon of Brownian motion Norbert Wiener [3j 
introduced as early as 1923 a certain probability distribution /j, concentrated 
on the set W d of continuous paths w : [0, oo[— * M. d , 1 1— > w(t) from the positive 
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half-line [0, oo[c M 1 into ci-dimensional Euclidean space M. d , d 6 {1, 2, 3, . . . }, 
which start at the origin, w(0) = 0. This distribution is a mathematically 
well-defined positive measure in the sense of general measure theory [3] and 
therefore induces a corresponding concept of integration over paths, which 
we denote by f wd fi(dw)(-). The (standard) Wiener measure (jl is uniquely 
determined by requiring that it is Gaussian with normalization 



for all j,k £ {l,...,d} and all t,s G [0, oo[. Here Wj(t) denotes the j- 
th component of the path w evaluated at (time) parameter t > 0. The 
simple Wiener integrals (QjJ and (|2|) imply that the components of w are, in 
probabilistic language, centered, independent and identically distributed. 

To Wiener integrals apply, in contrast to Feynman path integrals, all the 
rules and computational tools provided by general measure and integration 
theory, most notably Lebesgue's dominated-convergence theorem [4|. As a 
consequence, Wiener integration often serves, via the Feynman-Kac formula, 
as an efficient technique for obtaining results in quantum mechanics with 
complete rigour. An impressive compilation of such results was given by 
Barry Simon already in 1979. Since then not much has changed which is 
reflected by the fact that the second edition of his book [5] differs from 
the first one only by an addition of bibliographic notes on some of the 
more recent developments. Still Wiener integration should be considered 
neither as a secret weapon nor as a panacea for obtaining rigorous results 
in quantum mechanics. In any case, the Feynman-Kac formula is more 
than just a poetic rewriting of a Lie-Trotter formula. Ironically, Richard 
Feynman himself took advantage of that as early as 1955 in his celebrated 
paper on the polaron [6], in particular, by using the Jensen inequality. 

Now, what is the Feynman-Kac formula? Let us consider a spinless 
charged particle with configuration space M. d subjected to a scalar potential 
v : R d — > M},q I— ► v(q),q = {qi, ■ ■ ■ , qd) and a vector potential a : R d — > 
R , q i— > a(q),a = (ax, . . . , a^). The latter generates a magnetic field (tensor) 
defined by bjk := da^jdqj — daj/dq^- The corresponding (non-relativistic) 
quantum system is informally given by the Hamiltonian 




(1) 



and first, respectively, second moments given by 





(3) 
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where Q = (Qi, ■ ■ ■ ,Qd) and P = (Pi,..., Pa) denote the d-component 
operators of position and canonical momentum, respectively. They obey the 
canonical commutation relations QjPk—PkQj = iMjfcl. Here i = y/—l is the 
imaginary unit and h > is Planck's constant (divided by 2ir). Moreover, we 
have chosen physical units where both the mass and the charge of the particle 
are equal to 1. Under rather weak assumptions [7J on a and v, H(a,v) can 
be defined as a self-adjoint operator on the Hilbert space L 2 (R d ) of all (equi- 
valence classes of) Lebesgue square- integr able complex- valued functions on 
M d . Furthermore its "Boltzmann-Gibbs operator" e~P H ( a > v ) even posesses 
for each (3 e]0, oo[ an integral kernel (q\e~P H ( a,v ' \q') (in other words, position 
representation or Euclidian propagator) which is jointly continuous in q, q' £ 
R d . 

We are now prepared to state the Feynman-Kac formula. Apart from 
mathematical subtleties its content is most concisely expressed by the fol- 
lowing representation free version: 



e -PH(a,v) = I ^ dw ) e -MP^)-P/h e ^ p ( _ f dri;(u;(rn 2 )l + Q)} 
Jw d ^ Jo > 

xexpji/i J dTw(Th 2 ) • a(w(Th 2 )l + Q)} 



(4) 



Here the dot " • " between two d-component quantities refers to the Eu- 
clidian scalar product of R d and the integral containing the vector poten- 
tial is a suggestive notation for a stochastic line integral in the sense of 
R. L. Stratonovich and D. L. Fisk (corresponding to a mid-point discretiza- 
tion) . In the Wiener integration serves to disentangle the non-commuting 
operators P and Q in e~^ H ( a ' v ). For related remarks see Ref. [8] and ref- 
erences therein. We recall from (J2j> that fi neither depends on (3 nor on h. 
Also P/H is independent of H. 

By going informally to the position representation of ([!]) one gets the 
rigorously proven formula 

{q\e-P HM \q') = / n(dw)d(w( y (3h 2 )+q' -q) exp { - / dTv(w(rH 2 )+iA} 
Jw d L Jo J 

xexp^ihj dTw(Th 2 ) • a(w(Th 2 ) + q')Y (5) 

It even holds for a class of potentials v for which H (a, v) is not bounded 
from below [7] and remains true (by the self-adjointness of H (a, v)), if on its 
right side q and q' are exchanged and simultaneously i is changed to — i. The 
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Dirac delta in ([5]) indicates that all paths w to be integrated over arrive in 
q-q' £R d at "time" (3h 2 > 0. In fact, they may be considered to end there, 
because /i is Markovian and the Wiener integrand in ([5]) does not depend on 
w(tH 2 ) for t > j3. More precisely, the path integration may be performed 
with respect to the Brownian bridge [51 E]. In the Appendix below we shall 
present what we think is an illuminating derivation of the "bridge version" 
of ©, although in the following we shall not make (explicit) use of that 
version. 

In the next three sections we are going to illustrate the usefulness of ([5]) 
by deriving three inequalities of interest in quantum mechanics. 



2 Diamagnetic inequality 



Theorem 2.1. 

(g|e-^°' w V) 
holds for all P > and all q,q' 6R° 



< <9|e 



-PH(0,v)\ n i\ 



(6) 



Proof. Inequality Q is an immediate consequence of ([5]) by taking the ab- 
solute value, applying the "triangle inequality" 



(j,(dw) 



< I Kdw) |(-)| 
'w d 



and using the elementary identity |e s+1J '| = e x for x,y 6 M 1 . 
Remarks: 



(7) 
□ 



(i) This elegant proof is due to Edward Nelson, see Ref. [H] (also for other 
historical aspects of ([6]) and related inequalities). 

(ii) If J Rd dqe~ l3v ^ < do, the free energy -f3~ l In J Rd dg {q\e~^ H ^' v "> \q) at 
inverse temperature (3 > exists and ([6]) then implies that it cannot 
be lowered by turning on a magnetic field. Under weaker assumptions 
on v, for example for the hydrogen atom (that is, for d = 3 and v{q) = 

with 7 > 0) (0) still implies in the limit (3 — > oo the same sort 
of stability for the ground-state energy. Altogether this explains the 
name diamagnetic inequality. 

(iii) There are also diamagnetic inequalities in case the particle is restricted 
to a region in M. d of finite volume with Dirichlet, Neumann or other 
boundary conditions \10\ lllj. Moreover, the proof of the diamagnetic 
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inequality easily extends to the case of many (interacting) particles, 
provided there is no spin and no Fermi statistics involved. 

An interesting question is what can be said if a ^ is changed (pointwise) 
to another vector potential a' 7^ 0. For a partial answer see Sec. H] below. 

3 Quasi-classical upper bound on the integrated 
density of states in the case of a random scalar 
potential 

In the single-particle theory of electronic properties of disordered or amor- 
phous solids the scalar potential v in H(a, v) is considered to be a realization 
of a random field on M. d which is distributed according to some probabil- 
ity measure v on some set V of potentials v. We denote by j v v{dv) (•) 
the corresponding (functional) integration or averaging. One example is 
a Gaussian v with vanishing first moments and second moments given by 
fyv(dv) v(q)v(q f ) = C(q — q') for all q, q' S R d with some (even) covariance 
function C : R d — > R 1 . The fact that the second moments only depend on 
the difference q — q' reflects the assumed "homogeneity on average" . We also 
assume that C is continuous, C(q) tends to zero as \q\ -^00 and the single- 
site variance obeys < C(0) < 00. The R rf - homogeneity together with the 
decay of the correlations of the fluctuations at different sites with increasing 
distance implies the M. d -ergodicity of the (Gaussian) random potential. 

A quantity of basic interest in the above-mentioned theory is the inte- 
grated density of states. It may be defined |T2l [131 [7] as the non-decreasing 
function N : R 1 -> R 1 , E i-> N(E, a, q) where 

N(E, a, q) := [ v{dv) (q\Q(E - H(a, v))\q). (8) 

Here O denotes Heaviside's unit-step function and the (non-random) vector 
potential a as well as the position q G R rf are considered as parameters. If 
the random potential (characterized by v) and the magnetic field (generated 
by a) are both homogeneous, then N(E, a, q) actually does not depend on q. 
Of course, in the physically most relevant cases the random potential should 
be even ergodic, so that N(E, a, 0) coincides for ^-almost all realizations v 
with the number of eigenvalues per volume of a finite- volume restriction of 
H (a, v) below the energy E £ R 1 in the infinite-volume limit. Nevertheless, 
the following estimate holds also for random potentials and magnetic fields 
which are not homogeneous. 
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Theorem 3.1. If the probability measure v of the random potential has the 
propert 1 ^ that Lp := esssup rgIR d j y v{dv) e~^ v ^ < oo for all [3 > 0, then 

N(E, a, q) < (IvffiT^Lp ^ E (9) 
holds for all energies EgM 1 and all [3 > 0. 
Proof. 

N(E,a,q)e- (3E < [ u(dv) (q\e-^ H{a ' v) \q) < [ u(dv) (q\e- pHi °' v) \q) (10) 
Jv Jv 

Ii{dw)5{w{l3h 2 )) [ ^(dt;)exp{- f dr v(w(Th 2 ) + q)\ (11) 
rd Jv ^ Jo J 



!% [ ^dw)6(w((3h 2 )) ! u(d 
Jo P Jw d Jv 



< I I fi(dw)6(w((3h z )) I v {dv)e-^ w{Th2)+q) (12) 
Jw d 

< (2Til3h 2 y d/2 Lp. (13) 

Here (jlOp is due to the elementary inequality @(E — H(a,v)) < e@( E ~ H ( a > v >) , 
referring to the spectral theorem, and ([6]). Eq. © then gives (jlip . The next 
inequality is Jensen's with respect to the uniform average J^dr (•). The 
claim now follows from the definition of Lg, Eq. ([5]) with (a,v) = (0,0) and 
f£dr = p. The various interchanges of integrations can be justified by the 
Fubini-Tonelli theorem [4]. □ 

Remarks: 

(i) Theorem 13. 11 is a slight extension of a result which goes back to Pastur, 
see Thm. 9.1 in Ref. [12J. The right side of ([9]) is quasi- classical in the 
sense that it does not depend on a and does not take into account, 
due to the Jensen inequality (|13p in its proof, the non-commutativity 
of the kinetic and potential energy. 

(ii) While the estimate (J9|) holds for rather general random potentials, the 
various inequalities in its proof are responsible for its roughness, even 
when optimized with respect to (3 > 0. Nevertheless, it shows that 
N(E) decreases to at least exponentially fast as E — > —oo. For 
a homogeneous Gaussian random potential (and a constant magnetic 
field) the optimized estimate even reflects the exact Gaussian decay 



1 Given a function / : R d — > R 1 , then esssup rgR d |/(f)l denotes the smallest M £ [0, oo] 
such that |/(r)| < M holds for Lebesgue-almost all r € R d . 
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[121113] kiN(E) E 2 /2C{0) as E -> -oo. For a non-Gaussian 

random potential, like a repulsive Poissonian one, the leading low- 
energy decay of N, the Lifshitz tail, typically is of true quantum nature 
[T2| [13], [H] and can therefore not be reflected by the right side of 
Q. Although the (universal) leading high-energy growth N(E) ~ 
{E/2-Kh 2 ) d / 2 /r(l + d/2) as £ -> oo, see Refs. [121 [13] and references 
therein, is quasi-classical, the optimized right side of ([9|) overestimates 
it slightly by a constant factor (due to the elementary inequality in 

CUD). 

(hi) For related quasi- and pseudo-classical bounds (with a non-random v) 
we refer to Sec. 9 in Simon's book [5] and to Ref. [T5] . 



4 A simple diamagnetic monotonicity 

For a partial answer to the question raised at the end of Sec. [2] we only 
consider the planar case d = 2 with v = and a perpendicular magnetic 
field not depending on the second co-ordinate q 2 . We assume that b :=b\2 
is a continuously differentiable function of q± € Mr. One possible vector 
potential generating 6, not depending on q 2 either, is given by a^ b \q) := 
(0, Jq 1 dr 6(r)) . In the following theorem H^ b ' denotes any Hamiltonian on 

L 2 (M 2 ) which is gauge equivalent to H(a^,0) = if/2 + (P 2 - aJ, b) (Qi)) 2 /2 
for the given 6. The assertion (|14h is therefore gauge invariant. Nevertheless, 
in the proof we will use H(a^ b \0) and see that one can dispense with the 
absolute value on the right side of (I14h in this particular gauge. 

Theorem 4.1. If b and B are two magnetic fields as just described and 
satisfy either \b{r)\ < B(r) or \b(r)\ < —B{r) for all r £ Mr, then 

(q\e-P HiB) \q')\ < | ((q h 0)\e^ H{b) \(q[, 0))| e -(^) 2 /(^ 2 ) (14) 

ZioWs /or a// (3 > 0, all q = {q\, q 2 ) € M 2 and all q' = (q'i,q' 2 ) £ K 2 - 

Proof. By ([5]) the left side of (|14p is invariant under a global sign change 
of B. Therefore it suffices to consider the case B(r) > 0. For notational 
transparency we put h = 1 and write a and A for a 2 and j respectively. 
For a given pair ((3,w) G]0, oofxVF 1 we introduce the notations 

mp(a,w) := P 1 / (1to(w(t)), (15) 





4(a,«;) -/T 1 / dr(aKr))) 2 - (m /3 (a,^)) 2 (16) 

JO 
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for the mean and variance of a(u;(r)) with respect to the uniform aver- 
age Jo dr (•), and similarly with A instead of a. Next we observe the 
following two "doubling identities" 

2/3 2 [4(A,w)-4{a,w)] (17) 
P rP 

dr / dal [A(u;(t)) - A(w(<t))] 2 - [o(w(t)) - a(w(a))] 2 j (18) 

/3 rP 

dr da[a + (w(T)) - a + (w(a))][a^(w(T)) - a-(w(a))]. (19) 



The last integrand is non-negative, because the two functions r i— * a±(r) := 
A(r) db a(r) = Jq dr'(B(r') ± 6(r')), r £ R , are both non-decreasing since 
B(r') > \b(r')\ > =F6(r') by assumption. The same arguments apply when 
the path w is replaced by the rigidly shifted one w + q[ £ W 1 + M 1 defined 
by (w + 9i)(t) := w(r) + q^. To summarize, we have shown so far that 

s}(a,w + q[) < S }(A,w + fi). (20) 
Since H^ B \ in the particular gauge chosen, commutes with P2, we have 

= [ dkH( B \k)®\k)(k\, (21) 



using an informal notation for a direct-integral decomposition. Here the 
one-parameter family of effective Hamiltonians 

H^(k) := P 2 /2 + (kl - ,4(Qi)) 2 /2, kR 1 , (22) 

acts on the Hilbert space L 2 (M 1 ) of (wave) functions of the first co-ordinate. 
By ([21]) and ([22]) we get 



(q\e- f3H{B) \q'} = (27r)- 1 / dk^e'^^) (23) 
= (2tt/?)- 1 /V ((?2 -' ? 2 )2/(2/3) f fi(dw)8(w(fi) +q[- qi) 



JW 1 (24) 
x expj - /3sp(A,w + q[)/2\ exp|i(g 2 - q' 2 )mp(A, w + q[)\. 

Here we have used ([5]) with d = 1, a = and v = (k — A) 2 /2 and then 
performed the (Gaussian) integration with respect to fc. By applying the 
"triangle inequality" to (|24h and then using (1200 we finally obtain 
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,(l2-q' 2 ) 2 /(2P) 



(q\e p \q) 



(25) 



< (27T/3)- 1 / 2 j ^(dw)S(w(J3) +q[- qi) exp { - 0s%a, w + q[)/2] 

= (( qi ,0)\^ HQ,) UM = |((gi,0)|e-^ (6) |(^,0))| . (26) 
The last two equalities follow again from (124p with a instead of A. □ 
Remarks: 

(i) To our knowledge, Theorem 14.11 first appeared in Ref. |16j . It comple- 
ments some of the results obtained by Loss, Thaller and Erdos [171118]. 
For a survey of results of this genre see Sec. 9 in Ref. [19] . 



(ii) For a given sign-definite B the right side of (|14h can be made explicit 
by choosing for b the globally constant field Bq := inf rgK i |-B(r)|, so 
that 

q'^Bo/Ah (q 2 -q' 2 ) 2 



-fiH {B) I 



< 



Bq/AttH 



exp 



}■ 



smh((3hB /2) I tanh((3hB /2) 2(3h 2 

If Bq ^ 0, the Gaussian decay on the right side is faster along the 1- 
than along the 2-direction. Such an anisotropy has been found also for 
the almost-sure transport properties in the case that B is a (Gaussian) 
random field with non-zero mean |20| . 



Appendix 

For convenience of the reader we are going to derive the bridge version of 
the Feynman-Kac formula ([5]). We start out from a fixed triple (T,q,q') £ 
]0, oo[ xl^ x M. d and associate to each continuous path w : [0, oof— > W 1 with 
w{0) = another path w : [0, T] -► R d defined by 

w(t) :=w(t)+q'-^(w(T)+q'-q), t € [0,T\. (27) 

Obviously, w is a bridge path in the sense that w £ q , where 

n T,q, q > : = : [0, T) -> R d | w is continuous, w(0) = q',v(T) = g} (28) 

is the set of all continuous paths connecting position q' to position q in the 
time period T. In fact, Clj, , is the image of W d under the mapping w \—* w 
given by ([271) . 



9 



Writing ((•)) := J wd /jt(dw)(-) for the expectation or averaging induced by 
the Wiener measure, Eqs. ([1]) and ([2]) yield 

(® j (t))=</ j -~tf j -q j ), (29) 

{wj(t)w k (s)) - {wj(t))(w k (s)) = 5 jk (^mm{t,s} - (30) 

for the first and second moments of the bridge paths (|27|) . in other words, 
for their mean and covariance. 

Hence, by the (affine) linearity and the surjectivity of the mapping 
w i— » w the induced image of the Wiener measure /i on W d is a Gaus- 
sian probability measure pT,q,q' ° n , with mean and covariance given 
by the right sides of (f29|) and ([30]) . Other useful consequences of ([2]) and 
(j27|) are the two equalities 



(w fe (T)%(i)) = = (w k (T))(wj(t)). (31) 

They hold for all j,k E {1, ...,d} and all t E [0, T] and imply, by the 
Gaussian nature of the Wiener measure /U, that the random point w(T) E M d 
is independent of the family of random points {^(i)}te[o,T] C 

Now, how can all this be applied to the right side of ([5])? First we get 
from (2ZD that 

J(«;(T) + g 7 - g)F T (u» + g') = 5(u»(T) + g' - g)F r (w) (32) 

for any complex-valued function(al) Ft on W d + R rf , not depending on the 
points w(t) of w for t > T. The two mentioned independencies then give 

(5(w(T) +q'- q)F T (w + q')) = (5(w(T) + q' - q))(F T (w)) (33) 

= {2irT)- d / 2 e-(i-^ 2 / 2T (F T (w)) (34) 

and the bridge version of the Feynman-Kac formula © eventually reads 

( q \ e -P H M\q>) = (2TTf3h 2 )- d / 2 e-^ 2 / 2 ^ 2 



Qd 



(du>) exp 

x exp ^ih dr dj{rh 2 ) • a{uj{Th 2 )) }. (35) 



Remaining matters of rigour can be supplied. 

We note that a "one-parameter decomposition" similar to (|27p for a 
Gaussian random potential with non-negative covariance function has turned 
out to be useful in proving local Lipschitz continuity of the corresponding 
integrated density of states ([S]), see the proof of Cor. 4.3 in Ref. [TU] . 
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